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1 . 



Introducti on 



The study of atmospheric tidal oscillations has been the 
subject of long and intense research by many eminent scientists. 
Its history and theory have been particularly veil recounted by 
Chapman DO and ’.Vi Ikes DG • The regular 24-liourly fluctuations 
of pressure have been observed and analyzed since the invention 
of the barometer. More recently the air motions of the ionosphere 
have received special attention. Although solar diurnal variations 
of the winds in the troposphere and stratosphere have been the 
subject of inference from somewhat over-simplified theories, 
observations of such effects have been meager. 

The purpose of this study is to make an analysis of a 
possible solar semi-diurnal variation in the winds at ten 
kilometers from the regular upper wind observations published 
in the Northern Hemisphere Data Tabulations no . 
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Data selection 



The ten kilometer winds during a 100-day period from 
1 November 1959 through 8 February I960 observed at 15 stations 
were extracted from the Northern Hemisphere Data Tabulations DQ. 
Nine of these stations recorded the ten kilometer wind at 0000Z, 
0600Z, 1200Z, and 1800Z with reasonable regularity each day. The 
other six stations observed the ten kilometer wind at 0000Z and 
1200Z only. As a maximum of only four observations per day 
presents considerable impediment to satisfactory analysis of a 
semi-diurnal oscillation, the stations were selected in a spatial 
relationship designed to indicate the wave length of the oscilla- 
tion. Also a period was chosen to provide relative constancy of 
certain solar tidal parameters. The period selected is symmetrical 
to the northern hemisphere winter - the sun 1 s declination pro- 
gressing from about 15 degrees south on 1 November through the 
solstice on 22 December and back to about 15 degrees south on 
8 February. The location of the stations was such that they 
extended over approximately 180 degrees of longitude with the 
goal of encompassing spatially one complete semi-diurnal wave. 

As there is empirical evidence E3I1 f or a meridionally-directed 
standing pressure wave of solar semi-diurnal character with a 
node at 35 degrees latitude, the stations vrere selected primarily 
in the latitude belt 30 to 35 degrees north in order to minimize 
the wind effects of the standing wave. 
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3. Data analysis. 

The -westerly component of each six-hourly wind observation 
was c omputed' and* tabulated by station, date, and time. For each 
six-hourly time and at each station the average was taken. These 
results are tabulated in Table 1 along with the total number of 
observations comprising each average. 

Data from each of the nine stations recording four wind 
observations per day was treated first by the method of differ- 
ences between each pair of successive six-hourly averages. 
Assuming then that the change in westerly component over a 12- 
hour period is characteristic of a 24-hour wave and that one- 
half of this change occurs in each of the included six-hour 
periods, the diurnal effect is approximately eliminated by sub- 
tracting one-half of the 12 hour change from each of the six- 
hourly changes. This process' is applied over the periods 0000Z 
to 1200Z and from 1200Z to 0000Z. This method of analysis gives 
an apparent six-hourly change at each station of the form: 

00Z to 06Z : /^u = x ; 12Z to 18Z : ^ ■ u = y 

06Z to 12Z : AgU = -x ; 18Z to 00Z : u = -y 

Since x and y are in general unequal, the differences may be 
accounted for by a terdiurnal wave. For simplicity let this 
wave be symbolized by T, and assume it to be of constant ampli- 
tude and phase at each station during the observation period. 
Then the four differences just formed, +_ x and y, include 

o 

terdiurnal differences T^- T QQ , T^- T O0 , T lg - T J2 , and 
T 00” T 18 ‘ IIence Aj u i Ao u 5 A3 11 > 2111(3 A4 11 ma y be written in 
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the forms 

AjU = X' + (T 06 - T 00 ) 5 A 3 U = x' + (T 18 - T 12 ) 

Aj u = “’^ T ia“ T 08 O ’ 4i u " “ ( T 00 - T !bO 

where x' is the desired six-hour difference based on a semi- 
diurnal wave al one . Then may be formed 

V • A> u + V - A, u " 4 *' + 2<T 06 + T 18> " 2 (T 00 + T 12 ) (1> 

Since points on' the terdiurnal wave which are 12 hours apart are 

ISO degrees out of phase, it follows that T^q+ ^12 = ^ anc * 

T 10 = 0. Hence equation (l) eliminates the terdiurnal wave. 

Uo io 

Since the 24 hour wave had been previously eliminated, equation 
(l) gives four times the "amplitude" of the semi-diurnal change 
for the particular phase of the wave at the given station. In • 
other words, the mean magnitude of the six-hour difference at 
any station is 

A!" - A 2 “ + A 3 u - 4,“ 

4 *’ (2) 

and the four differences freed of both 24-hour and eight-hour 
harmonics are 

00Z to 06Z : x* ; 12Z to 18Z : x 1 

06Z to 12Z : -x 1 ; 1SZ to 00Z : -x ! 

The four differences shown above will be termed ^the "filtered" 
semi-diurnal six-hour differences. 

Oscillations of shorter periods which are evenly divisible 
into six have been removed also by the technique just described. 
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A four-hour period wave may not have been removed by this process; 
however, a four-hour wave of significant amplitude is unknown at 
the present tine. 

As the 100-day period selected is of such duration to com- 
prise three and one-half oscillations through the moon's range of 
declination and phase, lunar tidal effects are deemed to be small. 
For more nearly complete elimination of the lunar effect, selection 
of a period comprising at least 12 complete lunar oscillations 
( approximately 336 days ) would have been more appropriate; 
however, the lunar tidal oscillation is generally recognized to be 
small compared to the solar semi-diurnal wave. 

Appendix 1 contains the filtering calculations for each of 
the nine stations which, report four times daily. Figures 1 through 
9 depict graphically these filtered six-hour changes and the 
corresponding observed six-hourly means of the westerly wind 
component at each station. The filtered six-hour changes have been 
plotted in such a way as to divide the quantities x 5 , -x* 
symmetrically about the 100-day mean. Thus if a positive change 
occurs between 0000Z and 0600Z, a value x'/2 above the mean was 
plotted at 0600Z while at 0000Z a value -x’/2 was plotted. 
Examination of Figures 1 through 9 shows that the observed mean 
six-hour changes frequently have two maxima and two minima in the 
24-hour period and are in this respect and in essential form 
similar to the filtered six-hour differences unless the latter 
are so small as to be obscured by "noise 11 . 

A further analysis has been performed by plotting the 
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filtered fluctuations over t'ae full range of longitude of all 
nine stations. Figures 10 and 11 depict these spatially-arranged 
fluctuations for the periods 0000Z to 0600Z and 0600Z to 1200Z 

4 

respectively. For contrast, Figure 12 was prepared presenting 
the observed 12-hour changes for all 15 stations also plotted 
over the full, range of longitude. The curves inserted in these 
figures have been fitted by eye, and in all cases the vertical 
scale is the same to help convey relative magnitudes. 
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Fig. 1. Average westerly component of the wind at ten kilometers 
at Ship 034. 




Fig. 2. Average westerly component of the wind at ten kilometers 
at Norfolk 
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F ip;. 3. Average westerly component of the wind at ten kilometers 
a t l [ on tg on e r y ♦ 




Fig, 4, Average westerly component of the wind at ten kilometers 
at Shreveport. 



9 






Fig. 5. Average westerly component of the wind at ten kilometers 




Fig. 6. Average westerly component of the wind at ten kilometers 
at Santa Monica. 
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Fig. 7. Average westerly component of the wind at ten kilometers 
at Ship 134. 




at Midway. 
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Fig. 9. Average westerly component of the wind at ten kilometers 
at Ship 236. 
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Fig. 10. Filtered six-hour variations in the average -westerly 
component of the wind at ten kilometers for the periods OOOOZ 
to 0600Z and 1200Z to 1800Z. 
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component of the T/ind at ten kilometers for the periods 0600Z 
to 1200Z and 1800Z to 0000Z. This curve is a mirror-image of 
that in Fig. 10. - 
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Fig. 12. Observed 12-hour variation in the average westerly 
• component of the wind at ten kilometers for the period 0000Z 
to 1200Z. 
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4. Discussion of results. 

If the mean values of an element are available at each 

solar hour of the day, derivation through harmonic analysis of 

* 

the 24-liour wave, as well as the semi-diurnal wave and some 
higher frequency harmonics, is relatively simple. This has been 
done for the surface pressure wave ( see, for example, Chapman 
DO an d V/i Ikes DO ) and for other elements. As a further example, 
Greenhow and Neufeld have isolated a semi-diurnal oscillation 

of the zonal wind between 80 and 100 kilometers by observations 
of the drift of meteor trails. More recently Mantis DO has 
investigated the existence of both diurnal and semi-diurnal waves 
in the meridional component of the wind at 29 kilometers. Hi s 
data was based on drifts of constant pressure balloons. On the 
other hand Pekeris Qf] has presented a more sophisticated theory 
in which he assumes a temperature structure similar to that known 
to occur in the atmosphere. Using essentially the Pekeris model, 
Wilkes lias shown that a phase change of 130 degrees occurs in the 
pressure wave in the neighborhood of 30 kilometers; however, at 
ionospheric levels the phase of the semi-diurnal pressure wave 
is very nearly again in phase with the observed surface semi- 
diurnal wave, as is approximately required by the dynamo theory 
of magnetic variations. According to the results obtained in this 
paper, the change in phase (by 180 degrees) appears to have 
occurred as low as ten kilometers, and the semi-diurnal wave in 
the zonal component has an amplitude of two meters per second at 
this same level over oceanic stations. 
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In the case of most upper-air meteorological data, however, 
one can expect no more than four synoptic reports per day. Hence * 
the usual harmonic analysis procedures cannot be performed. For 
this reason, the technique of the filtered six-hour ohanges was 
devised utilizing the requirement of additional verification by 
proper spatial-temporal phase relationsliipto support the postu- 
late of a 12-hour wave at ten kilometers. 

Assuming a semi-diurnal wave at each station, the times of 
maxima or minima can be inferred from the filtered six-hour 
changes only to within +3 hours of the filtered maxima and ^ 
minima. The best that can be said of the amplitude from a single 
station record is that it should be not less than half the 
filtered change. Although only a few of the nine stations appear 
to have significant filtered six-liour changes, the spatial con- 
figuration of these wind variations depicted in Figures 10 and 
H lends credence to their semi-diurnal character and extent. 

There appears to be a fairly good indication that the adjacent 
maxima and minima are very close to 90 degrees of longitude or 
six hours apart. Except for the observations at Midway, the 
amplitudes and phases over the oceans suggest a sinusoidal 
character. V possible explanation for the much smaller amplitude 
over the United States may be the presence of more intense 
vertical motions associated with persistent long wave troughs 
frequent in winter there, which could obscure the horizontal 
fluctuations sought. The results obtained at Midway in opposition 
to other indications may be due to its departure in latitude by 
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five degrees from the mean latitude of all the stations. This 



small latitude difference, however, does not strongly support 

such an explanation. 

* 

The phase of the wind fluctuation depleted in Figures 10 
and 11 deserves some comment because it is opposite to that 
described in the wcl 1 -known text by Mitra QQ • gi ves a 
simplified theory using an atmosphere of constant density. 
According to this theory the surface semi-diurnal wind is in 
phase with the semi-diurnal pressure wave, having a maximum 
easterly component at approximately ten a.m. and ten p.m. local 
time. This would mean that the prevailing westerlies at 33 
degrees north latitude should be diminished at these local times 
Figures 10 and 11, however, show the westerly wind to be greater 
than the mean during the period 0600 to 1200 local time. 

A statistical verification of certain of these results 
was considered in the following manner: Assuming a solar diurnal 
wave of 24-hour period of similar form as the assumed terdiurnal 
wave and symbolizing this wave with n S M , that portion of the 
westerly component of the wind due to tidal effects may be con- 
sidered the sura of x,., the semi-diurnal wave, T, the terdiurnal 
wave, and S. This may be stated as follows: 

U 00Z = X i,00Z + T 00Z + S 00Z + u > and 



U 12Z ' X i,12Z + T 12Z + S 12Z + U 



where u is the mean zonal wind for the period of study. Adding 
yields 
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U 00Z + U 12Z = X i,OOZ + X i,12Z + T 00Z + T 12Z + S 00Z + S 12Z + ?U ’ 

but T 00Z + T 12Z = ° aS bef ° re ’ S 00Z + S 12Z = ° alS0 ‘ The 
neon of the semi-diurnal components 12 hours apart thus corre- 
sponds to the mean of the zonal winds taken 12 hours apart. Thus 

•-( X i,00Z + X i,127.) = '^ U 00z + U 12Z^ “ U » 
and similarly 

“f X i,06Z + X i,18Z^ = “^ U 06Z + U 18Z^ " U * 

In other words the daily mean of the 0000Z and 1200Z zonal wind 

observations may be considered to be a point on the semi-diurnal 
.wave one-half period out of phase with the mean of the 0600Z and 
1800Z observations. The assumption must be made, however, that 
the number of 0000Z and 1200Z observations are equal (or nearly/ 
so), and that the same is true for the number of 0600Z and 1800Z 
reports. As the greatest semi-diurnal variation in the westerly 
component of the wind appeared to exist over the ocean stations, 
the Student*s t test was applied at Ships 034, / 134, and 236 to 
test the hypothesis that the means of the 0000Z and 1200Z obser- 
vations were significantly different from the means of the 0600Z 
and 1800Z observations as the filtering analysis had indicated. 
Table 1 shows that the number of observations at each of these 
stations reasonably fulfills the required condition of equality. 
Utilizing the symbolism and method of Iloel DO ■whe re 
x = the 0000Z and 1200Z observations, 
y = the 0600Z and 1S00Z obserCati ons , 
n = the number of observations, 
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II 



-2 

x = the square of the variance, and 



2 

s 



V 



t 



= n + n - 2 = 

x 7 . 
x “ 7 




2 

s" 

x 



2 

+ n s 

y y 



degrees of freedom, 

V n n (n +"n - 2?) 

JEJLJE I 1 

n 4- n 

x y 



Computation of t yields the following results with correspondin 

* 

probability that the means of x and y above are significantly 
different: 

Ship 034 : t = 1.15 : 87$ , 

Ship 134 : t = 2.16 : 98$ , and 

Ship 236 : t = 3.34 : 99$ . 

These results lead to the conclusion that the apparently signif 

icant difference between 'wind observations at 0000Z and 1200Z 
and those at 0600Z and 1800Z is real and may be due to a semi- 



diurnal' effect. 



5. Summary and conclusions. 

Previous studies calculate an expected semi-diurnal varia- 
tion of winds at the surface of the order of one-half meter per 
second. There appears* f rom this investigation suggesti ons of much 
greater semi-diurnal variations in the winds at the ordinary 
synoptic levels, particularly over the ocean areas. Such magni- 
tudes may require semi-diurnal modification of vertical shear 
prognostications so important in rocketry. This and other poss- 
ible applications may be reason for more extensive study of an 
apparently neglected meteorological parameter. 

Though accurate harmonic analysis of the various tidal 
effects requires more or less continous observation, considera- 
tion of the simultaneous combination of both spatial and temporal 
features of synoptical ly observed parameters may yield signifi- 
cant results. 
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APPENDIX l 



FILTERING CALCULA.T I ONS 



The calculations involved in the analysis of the observed 
six-hour changes in the -westerly wind component are completed 
as follows: 

Let u = average -westerly component of the wind. 

U 06Z ” U 00Z = Al U 5 U 18Z “ U 12Z = 43 U ’ 

U 12Z “ U 06Z = A 2 U 5 U 00Z “ U 18Z = * 

Elimination of the 24-hour wave is accomplished as follows; 



A u + A> u 
A U - ( 2 

Ai* + Ag u 
^2 U " ( 2 



) = 



) = 



x ; 
-x; 



Ai u 

Ai u 



(^-) - y , 



A" + ^4 U \ 

( - ) = _ y . 



Elimination of the eight-hour wave is accomplished as follows: 

Ai u - a 2 u + a 3 u - a* u 

This computation resulting in thevfiltered semi-diurnal 
wave may be plotted in station model form as follows: 



u ooz 


U 12Z 


U 06Z 


U 18Z 


A x u 






4j u 


Aj u + A 2 u 


Ai u + 4* u 


2 


2 


x . 


_JL_ 


-X 


'-y 


X' 


X 1 


-X' 


-X ' 
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Specific calculations for each of the nine stations 
recording the ten kilometer wind four times per day follow. 



Ship 


034 


Norfolk 


308 


Montgomery 


226 


35. ON 


48.0 V/ 


36. 9N 


76.2 V/ 


32. 3N 


86.4 V/ 


20 • o G 


18.69 


37.71 


37.05 


33.77 


33.04 


18.18 


14.66 


37.26 


38.61 


35.17 


33.66 


-2.35 


-4.03 


-0.45 


+1.56 


+1.40 


+0.62 


+0.51 


+5.87 


-0.21 


-0.90 


-2.13 


+0.11 


-0.92 


+0.92 


-O'. 33 


+0.33 


-0.36 


+0.36 


-1.43 


-4.95 


-0.12 


+1.23 


+1.76 


+0.26 


+1 . 43 


+4.95 


+0.12 


-1.23 


-1.77 


-0.25 


-3.14 


-3.14 


+0.56 


+0.56 


+1.01 


+1.01 


+3.14 


+3.14 


-0.56 


-0.56 


-1.01 


-1.01 



Shreveport 


24S 


Fort Worth 


259 


Santa 


Monica 295 


32. 5N 


93.817 


32 . SN 


97 . 117 


33. 9N 


118.417 


31.32 


31.41 


31.15 


31.36 


15.97 


18.73 


30.30 


30.51 


31.44 


29.72 


17.31. 


16.66 


-1.02 


-0.90 


+0.29 


-1.64 


+1.34 


-2.07 


+1.11 


+0 . 81 


-0.08 


+1.43- 


+1.42 


-0.69 


+0.05 


-0.05 


+0.11 


-0.11 


+1.38 


-1.38 


-1.07 


-0 o S5 


+0.18 


-1.53 


-0.04 


-0.69 


+1.06 


+0.86 


-0.19 


+1.54 


+0.04 


+0.69 


-0.96 


-0.96 


-0.63 


-0.63 


-0.37 


-0.37 


+0.96 


+0.96 


+0.63 


+0.63 


+0.37 


+0.37 
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Ship 


134 


Miclvay 


066 


Ship 


.236 


30. ON 


140.0V/ 


28. 2N 


177.4V/ 


34 . ON 


164. 0E 


14.41 


15.40 


22.27 


24.32 


53.80 


51.83 


10.94 


8.03 . 


24.78 


25.78 


40.27 


45.27 


-0.4? 


~7,a7 


+12.01 


+ 1.40 


-1.0. fi3 


—0 • 80 


+4.46 


' +6.38 


-0.46 


-3.51 


+11.56 


+8.23 


+0.50 


-0.50 


+ 1.03 


-1.03 


-0.99 


+0.99 


-3.97 


-C.87 


+1.48 


+2.49 


-12.54 


-7.25 


+3.96 


+6.8S 


-1.49 


-2.48 


+12.55 


+7.24 


-5.42 


-5.42 


+1.99 


+ 1.99 


-9.90 * 


-9.90 


+5.42 


+5.42 


-1.99 


-1.99 


+9.90 


+9.90 



) 



Additional stations recording the -wind at 0000Z and 
1200Z only are plotted in model as follows: 

U 00Z U 12Z 





12Z 


00Z 


OOZ 


12Z 




ICenitra 
34. 3N 


119 

6.6?/ 


Charleston 

32.-9N 


00 

O 

• ro 
o o 

-4 GO 


Tuscon 
32. IN 


274 

110.9V/ 


15.79 


15.72 


31.30 


31.91 


24.81 


24.60 


1 

o 

• 

o 

-5 


+0.07 


+0.61 


-0.61 


-0.21 


+0.21 


San Diego 
32. 7N 


290 

117.2V/ 


llil o 285 

19. 7N 


155. IV/ 


Wake 
19. 3N 


245 

166. 6E 


17.56 


20.61 


17.19 


17.62 


3.44 


3.71 


+ 3.05 


-3.05 


+0.43 


-0.43 


+0.27 


-0.27 
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